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Abstract. In this paper, we consider two regularized transformation-optics 
cloaking schemes for electromagnetic (EM) waves. Both schemes are based on 
the blowup construction with the generating sets being, respectively, a generic 
curve and a planar subset. We derive sharp asymptotic estimates in assessing 
the cloaking performances of the two constructions in terms of the regularization 
parameters and the geometries of the cloaking devices. The first construction 
yields an approximate full-cloak, whereas the second construction yields an ap¬ 
proximate partial-cloak. Moreover, by incorporating properly chosen conduct¬ 
ing layers, both cloaking constructions are capable of nearly cloaking arbitrary 
EM contents. This work complements the existing results in [5||^ on approx¬ 
imate EM cloaks with the generating set being a singular point, and it also 
extends on regularized full and partial cloaks for acoustic waves governed 

by the Helmholtz system to the more challenging EM case governed by the full 
Maxwell system. 
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1. Introduction 


1.1. Background. This paper concerns with invisibility cloaking of electromag¬ 
netic (EM) waves via the approach of transformation optics. This method which 
is based on the transformation properties of the optical material parameters and 
the invariance properties of the equations modeling the wave phenomena, was pi¬ 
oneered in 16 for electrostatics. For Maxwell’s system, the same transformation 
optics approach was developed in 34 . In two dimensions, Leonhardt 


24 used the 


conformal mapping to cloak rays in the geometric optics approximation. The ob¬ 
tained cloaking materials are called metamaterials. The metamaterials proposed 
in [T7l|24l|^ for the ideal cloaks are singular, and this has aroused great interest in 


the literature to deal with the singular structures. In 12,30 , the authors proposed 


to consider the hnite-energy solutions from singularly weighted Sobolev spaces for 
the underlying singular PDEs, and both acoustic cloaking and EM cloaking were 
treated. In 11, [^,20, 21,27 , the authors proposed to avoid the singular structures 
by incorporating regularization into the cloaking construction, and instead of ideal 
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invisibility cloak, one considers approximate/near invisibility cloak. The latter ap¬ 
proach has been further investigated in for the conductivity equation and 

Helmholtz system, modeling electric impedance tomography (EIT) and acoustic 
wave scattering, respectively; and in [5]-[^ for the Maxwell system, modeling the 
EM wave scattering. For all of the above mentioned work on regularized approx¬ 
imate cloaks, the generating set is a singular point, and one always achieves the 
full-cloak; that is, the invisibility is attainable for detecting waves coming from 
every possible incident/impinging direction, and observations made at every pos¬ 
sible angle. In a recent article [^, the authors proposed to study regularized 
partial/customized cloaks for acoustic waves; that is, the invisibility is only at¬ 
tainable for limited/customized apertures of incidence and observation. The key 
idea is to properly choose the generating set for the blowup construction of the 
cloaking device. In [^, the authors only proved qualitative convergence for the 
proposed partial/customized cloaking construction, and the corresponding result 
was further quantihed in by the authors of the current article. 

In this paper, we shall extend [9,25 on the regularized partial/customized cloaks 
for acoustic wave scattering to the case of electromagnetic wave scattering. We 
present two near-cloaking schemes with the generating sets being a generic curve 
or a planar subset, respectively. It is shown that the first scheme yields an ap¬ 
proximate full-cloak, whose invisibility effect is attainable for the whole aperture 
of incidence and observation angles. The result obtained complements |5]-[^ on 
approximate full cloaks. However, the generating set for the blowup construction 
considered in |5]-[^ is a singular point, whereas in this study, the generating set is a 
generic curve. The cloaking material in our full-cloaking scheme is less “singular” 
than those in [5]-[^ , but at the cost of losing some degree of accuracy on the invisi¬ 
bility approximation. The second scheme would yield an approximate partial-cloak 
with limited apertures of incidence and observation. We derive sharp asymptotic 
estimates in assessing the cloaking performances in terms of the regularization 
parameters and the geometries of the generating sets. The estimates are indepen¬ 
dent of the EM contents being cloaked, which means that the proposed cloaking 
schemes are capable of nearly cloaking arbitrary EM objects. Compared to |5||^ 
on approximate full-cloaks, we need to deal with anisotropic geometries, whereas 
compared to [9,25 on approximate partial-cloaks for acoustic scattering governed 
by the Helmholtz system, we need to tackle the more challenging Maxwell system. 
Finally, we refer the readers to [8,14,15,29,37 for surveys on the theoretical and 
experimental progress on transformation-optics cloaking in the literature. 


1.2. Mathematical formulation. Consider a homogeneous space with the (nor¬ 
malized) EM medium parameters described by the electric permittivity Sq = Isxs 
and magnetic permeability /iq = Isxs- Here and also in what follows, I 3 X 3 denotes 
the identity matrix in For notational convenience, we also let ao := 0 ■ Isxs 

denote the conductivity tensor of the homogeneous background space. We shall 
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consider the invisibility cloaking in the homogeneons space described above. Fol¬ 
lowing the spirt in [5]-[^, the proposed cloaking device is compactly snpported in 
a bonnded domain and takes a three-layered structnre. Let d fig ^ be 
bonnded domains snch that ^2^, f2c\^a and 12\12c are connected, and they repre¬ 
sent, respectively, the cloaked region, conducting layer and cloaking layer of the 
proposed cloaking device. Let Fq be a bounded open set in and it shall be 
referred to as a generating set in the following. For S G M+, we let Ds denote an 
open neighborhood of Fq such that Ds ^ Tq (in the sense of Hausdorff distance) as 
6 —?• -|-0. Ds will be referred to as the virtual domain, and shall be specihed below. 
Throughout, we assume that there exists a bi-Lipschitz and orientation-preserving 
mapping —)■ such that 


Fs{,Ds/ 2 ) = h2a, Fs{Ds\Ds/ 2 ) = ^c\^a, Fsi^XDs) = and = Identity. 

( 1 . 1 ) 

Next, we describe the EM medium parameter distributions cr} in the 

physical space containing the cloaking device, and cr^} in the virtual 

space containing the virtual domain. The EM medium parameters are all assumed 
to be symmetric-positive-dehnite-matrix valued functions, and they characterize, 
respectively, the electric permittivity, magnetic permeability and electric conduc¬ 
tivity. In what follows, e,/r, a} and /r^, cr^} shall be referred to as, 

respectively, the physical and virtual scattering conhgurations. Let 


and 


{R“;£, ftcr} 




hO) O'o 

in 



in 



in 



in 


^0, ho, O'o 

in 

M3 \ F,, 

Ph O'l 

in 

Ds \ Ds/2 

pa: O^a 

in 

DsI2- 


( 1 . 2 ) 


(1.3) 


The virtual and physical scattering conhgurations are connected by the so-called 
push-forward via the (blowup) transformation in (1.1). To that end, we next 
introduce the push-forward of EM mediums. Let m and ms, respectively, denote 
the physical and virtual parameter tensors, where m = e,p, or a. Dehne the 
push-forward {Fs)^ms as 


m = (Fs).ms := • (DFs f^ o Fp, (1.4) 

where DFs denotes the Jacobian matrix of the transformation Fs. Throughout 
the rest of our study, we assume that 

£, /r, a} = {Fs)^{R^-, es, p-s, os} ■= (F^)*^^, {Fs)^ps, {Fs)*as}. (1.5) 
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Next, we consider the time-harmonic EM wave scattering in the physical space. 
Let 

E‘(x) := HM=-(V X E*)(x), x G (1.6) 

UJ 

where p G ]R^\{0}, d G and u G M+. (E*,H*) in (1.6) is called a pair of 
EM plane waves with E* the electric held and H* the magnetic held, p is the 
polarization tensor, d is the incident direction and uj is the wavennmber of the 
plane waves E* and H*. There always holds that 

p T d, namely p ■ d = 0. (1-7) 

E* and H* are entire solntions to the following Maxwell eqnations 

V X E* — zca/ioH* = 0 in 

V X H* -|- = 0 in 

The EM scattering in the physical space e,/i, cr} dne to the incident plane 
waves (EbH*) is described by the following Maxwell system 

V X E — icn/iH = 0 in . 

V X H + i(n(e + i5)E = 0 in 

snbject to the Silver-Muller radiation condition: 

(H-H*) X X- (E-E*)) = 0, (1.9) 


lim 

llxll^oo 


where x = x/||x|| for x G ]R^\{0}. We seek solntions E, H G Hioc{cut\;R^) to 
(1.8); see for the well-posedness of the scattering system (|L^). Here 

and also in what follows, we shall often nse the spaces 

Hioc{cut\] X) = {U\b G iL(cnrl; B)] B is any bonnded snbdomain of X} 


and 


H{cm\]B) = {H G {L\B)Y- V x f/ G {L\B)Y]. 


It is known that the solntion E to (1.8) admits the following asymptotic expan¬ 
sion as ||x|| —)■ oo (see, e.g., [To]) 

1 


Efx) - ETx) = 


where 


Aoo(x;p,d) := 


; EM T O 


-;E* 


( 1 . 10 ) 


( 1 . 11 ) 


is known as the scattering amplitude and x denotes the direction of observation. 
It is readily verihed that 


Aoo(x; p, d) = IIpII Aoo ( x; -||^, d ) . 


( 1 . 12 ) 


Therefore, withont loss of generality and thronghont the rest of onr stndy, we shall 
assnme that ||p|| = 1, namely, p G 
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Definition 1.1. Let Sp C C and C 

near/approximate-cloak if 


{12;e,/i, cr} is said to be a 


|Aoo(x;p,d)|| <1 for x G E^, p G Ep, d G E^. 


(1.13) 


If Ep = Erf = E^ = then it is called an approximate fnll-cloak, otherwise it is 
called an approximate partial-cloak. 


According to Definitio n |l.l[ the cloaking layer {f2\f2c; /i*, cr/} together with 

the condncting layer {f2c\^a; } makes the target EM object {12^; ha; ^1} 

nearly invisible to detecting waves (1.6) with d G Erf and p G Ep, and observation 
in the apertnre E^. Erf and E^ shall be referred to as, respectively, the aper¬ 
tures of incidence and observation of the partial-cloaking device. For practical 
considerations, throughout the current study, we assume that the cloaking device 
is not object-dependent; that is, the cloaked content {Da) <^ 0 } is arbitrary 

but regular^ namely, are all arbitrary symmetric positive definite matri¬ 

ces. In Definition ( |1.13[ ) is rather qualitative, and in the subsequent study, 
we shall quantify the near-cloaking effect and derive sharp estimate in assessing 
the cloaking performance. To that end, the following theorem plays a critical role 
(cf. §0). 

Theorem 1.1. Let (E, H) G if;oc(curl; be the (unique) pair of solutions to 
(1.8). Define the pull-back fields by 

Erf = (Frf)*E := {DFsfE o Fs, Hrf = (Frf)*H := {DFsfU o Fs. 

Then the pull-back fields (Erf, Hrf) G iL^oc (curl; M^)^ satisfy the following Maxwell 
equations 


(1.14) 


V X Erf — iujfioHs = 0 

in 

M3\Drf 

V X Hrf -f icueoEs = 0 

in 

M3\Drf 

V X Erf — ica/XrfHrf = 0 

in 

Ds, 

V X Hrf -|- icn^^rf -\- i^^Erf = 0 

in 

Ds 


subject to the Silver-Muller radiation condition: 

(Hrf - H‘) X X - (Erf - E')) = 0, 


lim ||x 

llxll^OO 


Particularly, since Fs = Identity in ]R^\D, one has that 

A^i±:E,*) = i e 


(1.16) 


(1.16) 


where A^(x;E*) denotes the scattering amplitude corresponding to the Maxwell 
system dlTil ). 

Hence, by Theorem |1.1|, in order to assess the cloaking performance of the 


cloaking device {D; £,/i, a} in (1.2) associated with the physical scattering system 
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(1.8)-(1.9), it suffices for us to investigate the virtual scattering system (1.14)- 


(1.15). Here, it is noted for emphasis that by (1.5), one has 


{-D«/2;£a,;i.,CTa} = (U 


(1.17) 


Since {Ha; £*, /x*, cr*} is arbitrary and regular and Fs is bi-Lipschitz and orientation¬ 
preserving, it is clear that {Ds/ 2 ', £a, l^a, o'a} is also arbitrary and regular. 

In summarizing our discussion so far, in order to construct a near-cloaking de¬ 
vice, one needs to firstly select a suitable generating set Fq and the corresponding 


virtual domain Ds; and secondly the blowup transformation in (1.1); and thirdly 


the virtual conducting layer {Ds\Ds/ 2 ', fJ-h cri}', and hnally assess the cloaking 
performance by studying the virtual scattering system (1.14)-(1.15). In this ar¬ 


ticle, we shall be mainly concerned with extending the full- and partial-cloaking 
schemes proposed in [9,25 for acoustic waves to the more challenging case with 
electromagnetic waves. There the generating sets are either a generic curve or a 
planar subset, and the blowup transformations are constructed via a concatenating 
technique. Hence, in the current study, we shall mainly focus on properly design¬ 
ing the suitable conducting layers and then assessing the cloaking performances 


by studying the corresponding virtual scattering system (1.14)-(1.15). 


The rest of the paper is organized as follows. In Section 2, we collect some pre¬ 
liminary knowledge on boundary layer potentials, which shall be used throughout 
our study. Sections 3 and 4 are, respectively, devoted to the study of the full- and 
partial-cloaking schemes. 


2. Boundary layer potentials 

Our study shall heavily rely on the vectorial boundary integral operators for 
Maxwell’s equations. In this section, we review some of the important properties 
of the vectorial boundary integral operators for the later use. 

2.1. Definitions. Let D be a bounded domain in with a O^-smooth boundary 
dD and a connected complement M.^\D. Let Vod- denote the surface divergence 
on dD and H®(ciD) be the usual Sobolev space of order s G M on dD. Let u be the 
exterior unit normal vector to dD and denote by TH’^(ciD) := (a G z/-a = 

0}, the space of vectors tangential to dD which is a subset of W{dD)^. We also 
introduce the function space 

™^UdD) ; =|a G TR^dD)- Vazx ■ a G (^D)}, 


endowed with the norm 
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Next, we recall that, for oj G M+ U {0}, the fundamental outgoing solution 
to the PDO (A + in is given by 


a.hx) = 


47r||x| 


( 2 . 1 ) 


In what follows, if a; = 0 we simply write G^j as G. 

For a density function a G we define the vectorial single layer po¬ 

tential associated with the fundamental solution G^ introduced in (2.1) by 

^‘^[a](x) := [ G^{x-y)a{y)day, x G (2.2) 

JdD 

For a scalar density (p G the single layer potential is dehned similarly by 

5ob](x) := /" G^(x - y)(p(y) dcTj/, x G (2.3) 

JdD 

The following boundary operator shall also be needed 

lJr{dD) 


Al^[a](x) = p.v. z/x X V X / G^(x,y)a(y) day, 

JdD 


(2.4) 


where L‘^{dD) := TH°(cIZl), and p.v. signihes the Cauchy principle value. In 
what follows, we denote by Ad, <Sd and J^d the operators and 

respectively. 

2.2. Boundary integral identities. Here and throughout the rest of the paper, 
we make use of the following notation: for a function u defined on M.^\dD, we 
denote 


and 


du 

du 


uL(x) = lim m(x ± rz/(x)), x G dH, 

T- 4 - 1-0 


(x) = lim (Vxm(x ± Tuhc )), p(x)), x G dD, 

T^+O 


if the limits exist, where u is the unit outward normal vector to dD. 

It is known that the single layer potential satisfies the trace formula (cf. 


10,31 


= ±-/+(KS)' 


where (/C‘d)* is the L^-adjoint of /C^ and 


on dD, 


(2.5) 


IC^[a] := p.v. 


IdD 


d(Ca;(x-y) 
diy{y) 


ip{y) day, X G dD. 


The jump relations in the following proposition are also known (see 10,31 
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Proposition 2.1. Let a G Then .4.‘J^[a] is eontinuous on and its 

curl satisfies the following jump formula, 

1 / X V X .A‘^[a]|^ = =f| + Al‘^[a] ondD, (2.6) 

where 

i/(x) X V X |^(x) = lim z/(x) x V x .4,‘^[a](x ± ti/(x)), Vx G dD, 


Equipped with the above knowledge, the solution pair (E^, H^) in \ Di to 
(1.14) can be represented using the following integral ansatz, 

E 5 (x) = E*(x) + Vx X .A^Ja](x), X G \ Ds, 

(2.7) 

H 5 (x) = — Vx X E,(x) = H'(x) + — Vx X Vx X .A‘^Ja](x), x G 

ZUJ lUJ j' 

(2.8) 


_ 1 /o 

where by (2.6) the vectorial density function a G TH^V {dDfij satisfies 


2+-^Dj[a](x) = i^x(E,-E^)(x) 


, y G dDi. 


(2.9) 


3. Regularized full-cloaking of EM waves 

In this section, we consider a regularized full-cloaking scheme of the EM waves by 
taking the generating set to be a generic curve. As discussed in the introduction, 
this scheme was considered in our earlier work for acoustic waves. For self- 
containedness, we briefly discuss the generating set Eq and the virtual domain 
for the proposed cloaking scheme in the sequel, which can also be found in [^. Let 
Fq be a smooth simple and non-closed curve in with two endpoints, denoted 
by Pq and Qq, respectively. Denote by A^(x) the normal plane of the curve Fq 
at X G Fq. We note that N{Pq) and N{Qq) are, respectively, defined by the left 
and right limits along Fq. Let q G M+. For any x G Fq, we let ^g(x) denote the 
disk lying on -IV(x), centered at x and of radius q. It is assumed that there exists 
go G IR+ such that when g < go, ^q(x) intersects To only at x. We let be given 
as 

■= yq{x) X Fo(x), X G Fo, (3.1) 

where To is identihed with its parametric representation Fo(x); see Fig. for a 
schematic illustration. Clearly, the facade of denoted by and parallel to 
To, is given by 

:= {x g ■ n(x); X G Fq, n(x) G A^(x) fl §^}, (3.2) 

and the two end-surfaces of are the two disks ^q(Po) and ^g(Qo)- Let 
and be two simply connected sets with U S^qfiPo) and = 
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Figure 1. Schematic illustration of the domain Dq for the regular¬ 
ized full-cloak. 

Sq^ U ^qq{Qo), such that Sqq := U Sqg U 5'“^ is a C^-smooth boundary of the 
domain Bg^ :=£>“(,□ B^^ U For 0 < q < qo, we set 

- ^>0) + a = If- ■ (x - P„) + P„; X e b;| . 

Qo L 9^0 J 

and similarly, Bg := q/qo ■ {Bg^ — Qq) + Qq. Let Sg and Sg, respectively, denote the 
boundaries of Bg and Bg excluding and . Now, we set Bg := BgLiB^^UBg, 
and Sg := U SgU Sg = dBg. 

Henceforth, we we let 6 G M+ be the asymptotically small regularization param¬ 
eter and let Bs denote the virtual domain used for the blowup construction of the 
cloaking device. We also let 

Sj:=S/uSjuS(“. (3.3) 

denote the boundary surface of the virtual domain Bs- Without loss of generality, 
we assume that qo = 1. We shall drop the dependence on q if one takes q = 1. 
For example, B and S denote, respectively, Bg and Sg with q = 1. It is remarked 
that in all of our subsequent arguments, B can always be replaced by Bt-^ with 
0 < tq < qo being a hxed number. 

In what follows, if we utilize z to denote the space variable on Fq , then for 
every y G -D^, we define a new variable G Fq which is the projection of y 
onto Fq. Meanwhile, if y belongs to Bg (respectively Bg), then Zy is dehned to 
be Po (respectively Qo). Henceforth, we let ^ denote the arc-length parameter 
of Fq and 6, which ranges from 0 to 27r, be the angle of the point on N{^) with 
respect to the central point x(^) G Fq. Moreover, we assume that if 0 = 0, then 
the corresponding points are those lying on the line that connects ro(0 fo ri(^), 
where Fi is dehned to be 

Fi := {x -I- ni(x);x G ro,ni(x) G iV(x) fl — a hxed vector for a given x}. 
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With the above preparation, we introduce a blowup transformation which maps 
y E Ds to y E D as follows 


dl(y) = y:=^(y 


whereas 


dl(y) = y := 


y-Pp 

S 

y-Qo 

<5 


'y) + 

y e Dj, 

(3.4) 

+ Po, 

- + Qo, 

y e Di 
y e Dl 

(3.5) 


Next, we present the crucial design of the lossy layer in (1.3). Dehne the Jacobian 
matrix B by 

B(y) = VyJl(y), y e Ds. (3.6) 

Set the material parameters es, fJ-s and as in the lossy layer Ds \ Ds /2 to be 

^^(x) = £i(x) := /i5(x) =/ii(x) ^ J*|B|B"\ 

for xE Ds\ Ds/ 2 , 

stands for the determinant when 


cr5(x) = cri(x) := JbBlB 


-1 


where r, s and t are all real numbers and 
related to a square matrix. 

We are now in a position to present the main theorem on the approximate 
full-cloak constructed by using Ds described above as the virtual domain. 

Theorem 3.1. Let Ds be as describ ed abo ve with dDs = Ss defined in Let 

b e the pair of solutions to (l.lf), with {1 2; es , jis, C /i^, 

defined in (1.3), and {Ds\Ds/ 2 ',£s, h‘S,o's} given in (3.1). Define 

fi = min{l, —1 -I- r -I- s, —1 -I- f -I- s}, (3' = min{l, —2 -|- r -|- s, —2 -|- f J- s}. 

Ifr, s and t are chosen such that fi' — t/2 > 0, then there exists Jo ^ 1^+ such that 
when J < Jo; 

||Ai.(xip,d)|| <C(i'S-‘/2+i + i2) (3.8) 

where C is a positive constant depending on u and D, but independent of Sa, p-a, 
a a and x, p, d. 

Remark 3.1. Following our earlier discussion, one can immediately infer by Theo¬ 
rems 0 and |3.1| that the push-forwarded structure in ( |1.2[ ), 

{n;e,p,a} = {Fs)^{D.;es, ps, (^ 5 } 

produces an approximate full-cloaking device within at least J-accuracy to the 
ideal cloak. Indeed, if we set s = 2 and r = t = 0 then one has (3 = 1, fi' = 0 
and the accuracy of the ideal cloak will be which is the highest accuracy that 
one can obtain for such a construction. Particularly, it is emphasize that in (3.8), 
the estimate is independent of Pa, o'a, and this means that the cloaked content 
{12a; £*,/!*, cr*} in (1.2) can be arbitrary but regular. Finally, as remarked earlier, 
we refer to for the construction of the blowup transformation Fs, which we 
always assume the existence in the current study. 
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The subsequent three subsections are devoted to the proof of Theorem |3.1 
our later use, we hrst derive some critical lemmas. 


For 


3.1. Auxiliary lemmas. In this subsection we present some auxiliary lemmas 
that are essential for our analysis of the far held estimates. To begin with, we 
show the following properties of the blowup transformation dehned in (3.4) 


Lemma 3.1 (Lemma 4.1 in [^). Let A be the transformation introduced in (3.4) 
and (3.5) which maps the region Ds to D. Let B(y) be the corresponding Jacobian 
matrix of A{y) given by (3.6). Then we have 


B(y) = 


jl3x3 - (j - y e d {, 


Ij 

5-*-3x3) 


y e u Dh, 


(3.9) 


'S ^ 

where the superscript T denotes the transpose of a vector or a matrix. Furthermore, 


B(y)py = y e dDs, 

where Vy stands for the unit outward normal vector to dDs at y E dDs. 


(3.10) 


Remark 3.2. In view of the Jacobian matrix form (3.9), one can also hnd that the 
eigenvalues of B(x), x G are either 1 or l/J. Hence for any vector held V G 
there holds 

||Vf < (B(x)V,V) < J“iVf 


/ 


uniformly for x G D^. It can also be easily seen from (|3.9|) that 

|B(y)|=r^ yGl 
For the sake of simplicity, we dehne 

E + ■p TT'* ■tJ+ tJ Tji 

s .— rjs — iL, , Jrl^ .— Jrl^ — H. , 

Furthermore, we introduce the following notations 


m 




5- 


E(k) := E(A-i(k)) = E(x), H(k) ;= H(A-i(k)) = H(x), 
and dehne the corresponding helds after change of variables by 

E(i) := ((BL-‘E)(x), H(x) := ((BL'‘H)(x), 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


where 

B(i) := B(A-^(i)) = B(x). 

We mention that sometimes we write B and B in the sequel and omit their de¬ 
pendences for simplicity. The following lemma is of critical importance for our 
subsequent analysis. 
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Lemma 3.2 (Corollary 3.58 in (^). Let x = v4(x) with B(x) = VA(x). Then 
for the hounded domain Ds and any vector field V G H{cml; Ds), 

V(i) := (B^)-iV(i), V(i):=V(x), x G 

there hold the following identities 

|B|-iB(V X V){A-fi±)) = Vi X V(i), (3.16) 


and 




IdDs 


X V) 


Wda^ = 


>dD 


X V) ■ WrfcTj 


where W G H{cm\]Ds) and W{±) := (B^)-iW(i) := (B^)-iW(x). 
Note that 


day = 

with which one can show that 


S day, y £ s/, 
day, ye Sj U Sf, 


(3.17) 


(3.18) 


Lemma 3.3. Let D^, c G {/, a, b} be defined at the beg inni ng of this section. Let 
V and V, B be similarly defined as those in Lemma 
H{cutI;D), one has 


3.2 


Then for any W G 


'dDf 


'dD<^ 


i/i X V ■ W(5i;) da^ =6 
i/i X V ■ W(5i;) da^ =6 


'dDf 


ayyxY ■ W(x)da£, 


(3.19) 




i/jt X V ■ W(x) dax cG{a, 6}, (3.20) 


where W(i) := (B^)-^W(i) 


Proof. The proof follows directly from (3.17), (3.18) by using change of variables 
in the corresponding integrals. 


□ 


Lemma 3.4. Suppose a(x) = a(x) for x G dDs and x G dD. Define 


and 


Msf[a]{x) := p.v. z/i x Vz- x / G^(z£ - Zy)a(y) day, 

Jsi 


Al 5 c[a](x) := p.v. L-i X Vi X / G^(x - y)a(y) day. 


Define 

£,Dgi^) :=L'xXVxXVxX ^£)^[a](x). 

Then there hold the following results 

/ '5V<",[a|(x) + 0(32||al| TH-U2(aD))i ^ ^ 
^ \ Af5.[a](i)+ C>(5||a||TH-i/2(ai5)), x G 


(3.21) 
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and 


C 


DsWK^J 1 1,/- V Y7- V T7- V A -L n 


X G S'/, 

|Px X Vx X Vx X ^5.[a](x) + 0(5||a||TH-i/2(aD))) ^ ^ 


for c G {a, b}, where 


(3.22) 


/:g/[a](x) := Px X Vzj X Vzj; X / - Zji)a(y) da. 

Jsf 


yi 


and As<= is given in (2.2) by replacing dD and uj with S''^ and 0, respectively. 
Proof. For x, y G Ds, one has 

X - y = (x - Za,) - (y - Zj/) + Za, - Zy = 5 ((x - z^,) - (y - Zy)) + (z^, - Zy). 
Hence, we have the following expansion for x G S'/, 

||x - y|| = ||z£ - Zyll + 5((x - z^) - (y - Zy), ) + O(d^). 

II II 

Similarly 

(x - y, Px) = (zx - Zy, z/x) + 5((x - Zi) - (y - Zy), Vyf), 

l|x - y|r^ = ||zi - Zy\\-^ - 6 {{± - Zi) - (y - Zy), || J/_ ^Jia ) + 0 ( 6 ^), 


gi.;||x-y|| _ g*a;||zs-zg|| G ^ iujS{{± - Z^) - (y - Zy), 


Zx Zy 


\Zx Zy\ 


+ 0(5^). 


With those expansions at hand, we proceed to compute for x G S'/, 


VxGa;(x-y) = 


(x - y)e*"^ll^-yll 


47r||x-y||2 \^||x-y|| 


— lUJ 


47r||z£ - Zy||2 
= V.,G,,{zx-Zy) + 0{8), 
and by using vector calculus identities, 

Wl‘^Ja](x) = Px X Vx X / G^(x - y)a(y) day 

JdDs 


iuj + 0{5) 


= VxGa;(x - y)z/x ■ a(y) dcTy - / z/x ■ VxG'a;(x - y)a(y) dcxy 

JdDs JdDs 

=6 / V^-G^{zx - Zy)z/x ■ a(y) day - 6 z/x ■ V^-G^{zx - Zjj)a(y) day 
Jsf Jsf 

+ 0{6 ||a||rpjj-i/ 2 (-g^^) = ( 5 A 4 _g/[a](x) + 0{6 ||a||rpjj-i/ 2 (-g^)). 
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which proves the case of (3.21) for x G S'/. Next, note that if x G S*^, c G {a,b}, 
then 

Zi-Zy = 0 for y G S^, 

and 

II II All- -11 (S;-y,i^x) ^ 

l|x-y||=il|x-yl|. ||,_y|| - ■ x.yes,. 

With the above facts and by using a similar argument to the proof of the first case 
of (3.21), one can prove (3.21) for x G S|, c G {a,b}. To prove (3.22), we first note 
that 


V, X V, X ylS,[al(x) = x;"xlS,[al(x) + VW. • xl^,,|al(x). 
By using (4.17) in 1^, one can easily obtain 


(3.23) 




= ^ (3-24) 

J \ 5(^s- + ^^^)[a](x) + C>(52||a||TH-i/2(az))), x G S^", 

for c G {a, b}, where A‘^f is defined similarly as S‘^f in j^. Using similar expansion 
strategy, one can prove for x G S'/ 

VxVx ■ ^)(,Ja](x) = ■ ^^,[a](5c) + (3.25) 

Note that for x G S'^, c G {a, b}, there holds (see Proposition 5.2 in for details) 


VxVx ■ ^‘^Ja](x) = -VxVx ■ Asc[^m + (P(5||a| 


TH-U2(ao)J 


By combining (3.23)-(3.26) together, one can complete the proof. 


(3.26) 

□ 


3.2. Asymptotic expansions. In order to tackle the integral equation (2.9), 
we shall first derive some crucial asymptotic expansions. Henceforth, we denote 
a(y) := a(y) for y = A(y), y E Ds and y G dD. The same notation shall be 

adopted for x and x. _ 

For X G \ with sufficiently large ||x||, we can expand — E* from (2.7) 
in z G To as follows 


(E 5 -E‘)(x) = Vx A‘^Ja](x) = Vx [ 

JdD, 


y)a(y) day 


=V X 


'dDs 


G'a;(x - Zj^)a(y) day-V x 


dDs 

VG^(x - z 


I dDs 


y) 


(y - Zy)a(y) day 


+ V X ^ 
:=Ri + i?2 


(y-Zj/) V G^(x-C(y))(y-Zj^)a(y)day 


ODs 

+ R 3 , 


(3.27) 


where ({y) = rjy + {1 — rj)Zy G for some r] G (0,1), and the superscript T 
signifies the matrix transpose. We next estimate the three terms i?i, i ?2 and i ?3 in 
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(3.27). The term in (3.27) is a remainder term from the Taylor series expansion 


and it verihes the following estimate, 


<CS^ 


-R3||l°°(S2)3 — 

1 


Vx / (y-z^) V G^(x-C(y))(y-Zy)a(y)daj, 

JdDs 


|a||TH-i/2(g£))- 


L°°(S2)3 

(3.28) 


In the sequel, we shall need the expansion of the incident plane wave EMn z G Tq, 
and there holds 

OO 

E-(y) = E‘(z,) + VE‘(z ,). (y - z,) + - z»)”. (3.29) 

| a |=2 

where the multi-index a = {ai,a 2 ,a 3 ) and 9“ = dy^dy^^y^ with y = (i/i, 1/2, i/a)- 
Since for y G dDs, z/y = z/y, one further has 

OO 

z/y X E*(y) = UyxJ2 <5“a“E*(zy)(y - Zy)^. 

| o |=0 


In what follows, we dehne 

®(y) := ^(y) = ^yX E5(y) 


dDs 


(3.30) 


Theorem 3.2. Let E^ be the solution to {1.14h then there holds for x G 


/ ~ Zy)a{y) day =-2 ^^^(x - Zy)<I>(y) + C>((5(||<I)||TH-i/2(aD) + !))• 

Jsf Jsf 

(3.31) 

If one assumes that ‘h(y) = 0, y G dDs, then there holds for x G \ 71, 

(Es - E')(x) =2^2V X ( [ G,(x - zy)iyy x (VE'(z^)(y - Zy)) day 

^Jsf 

+ G^(x - z^)z/y X E\zy) day 

7s“us'> 

(VG^(x - Zy) ■ (y - Zy))vy X E*(zy)(y) da^ + 0{5^). 


isf 


Proof Recall that — | + -Mds is invertible on THjjy^(cR75) (see, e.g. 18 ). By 
(2.9), we see that 


a(x) = ( - 2 + -^Ds 


-1 


z/yX(E,-E*)(y) 


dDs-i 


X . 
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Using the results in Lemma 3.4 and the expansion of E* in (3.29), we have for 


X G S'/ that 


a(x) = - 2<h(x) + 2i/x X E*(zi.) - 45>l^^[$](x) + 2Sus, x (VE*(zi)(x - z,^)) 

+ 45A^g/[z/y X E (zy)](x) + (9((5 (||$||th-iG( 9£)) + !))• (3.32) 


Note that 


l•2^T 


I sf 


G^(x—Zy)(r'y xE*(zy)) dcTy = / G^{-x—Zy) / z/y (i-dxE*(zy) dr = 0, (3.33) 


'To 


where (r, i?) stand for the polar coordinates, and this together with (3.32) readily 
implies (|3.31l) . _ 


Next, if ‘h(y) = *h(y) = 0, then it can be seen from (3.21) and (3.32) that 


ITH 




where C* is a positive constant depending only on D and uj. By using our earlier 
results in (3.27) and (3.28), one can hrst show that 

(VGa;(x - Zy) ■ (y - Zj^))a(y) day + 0{5^) 
^^(x - Zy)a(y) day 


‘5Dja](x) 

=d 


G^(x - Zj^)a(y) day- 


dDs 


>dDs 


isf 


'sf 


G^(x - Zy)a.{y)day + 6" 

(VG^(x - Zy) ■ (y - Zy))a(y) day) + 0{S^). 


By using (3.32) we have for y G S'/ that 

a(y) = 2z/y X E'(zj^) + 2dz/y x (VE'(zj;)(y - Zy)) + 0(6^). 
Substituting (3.35) into ( |3.34[ ) and using (3.33) one can easily obtain 

‘5Dja](x) =26^ [ G^(x - z^)z/y x (VE*(zy)(y - Zy)) day 


(3.34) 

(3.35) 


Sf 


'S“US*> 


Guli'x. - Zy)z/y X E*(zjj) day 


'Sf 


(va 


y) ■ (y - Zy))i^y X E*(zy)(y) day) + 0{5^), 


which then completes the proof by using (2.7[). 


□ 


We continue with the estimates of Ri and R 2 in (3.27) for the proposed full¬ 
cloaking structure with arbitrary but regular fia and aa in Ds/ 2 - In what 
follows, we let G denote a generic positive constant. It may change from one 
inequality to another inequality in our estimates. Moreover, it may depend on 
different parameters, but it does not depend on pa, aa and p, d, x. 
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We first note that, by using (3.32) and (3.21) in Lemma 3.4, there holds 


ITH 


— ^(ll‘^llTH-U2(a£)) + !)• 


(3.36) 


Next, by taking expansion around z G Lq and using (3.31) one can show that for 
6 G M+ sufficiently small and ||x|| sufficiently large. 


l^ill = 
<5 


Vx / G'^(x - Zy)a(y) 

JdDs 


Vx / G'a;(x - z.y)a(y) dcTj; 


'sf 


+ 05“^ 


lTH-U2(a£)) 


<C6 


V X y ^ Gtj(x Zy)^{y)day +C5 ||^|| (||‘h||rpjj-i/2^g£,^ + 1) 


(3.37) 


—^^||x|| II*^IIth-i/2(5/) + C5'^ (||‘h||rpjj-i/2(g^) + 1). 


By using Taylor’s expansions again and (3.32), one can show the following estima¬ 
tion for i?2, 


\R2\\ = 


Vx / VG^(x-Zy) ■ (y - Zj/)a(y)daj 


'ODs 


<Cd^— , 11 , 1 


^IItH-i/2(s'/) + 1) + C'5 ll^ll (||4>||TH-i/2(g£)) + 1). (3. 


38) 


Hence, by applying the estimates in (3.28), (3.37) and (3.38) to (3.27), we have 
IIE5 — E II < C'n^||$||rpjj-i/2(5./^ -|- C6 j|^(||‘h||rpjj-i/2(g^) -|-1) (3.39) 


for ||x|| sufficiently large. With (3.39) at hand, one readily has 


Lemma 3.5. The scattering amplitude A^(x; p, d) corresponding to the scattering 


configuration described in Theorem 15*. j| satisfies 

ll"^oo(^j P) d)|| < C'(5||$||rpjj-i/2(g/) -|- (||^||TH-U2(-g^) + 1), 
where C depends only on D and oj. 


(3.40) 


3.3. Proof of Theorem |3.1[ By Lemma |3.5[ it is straightforward to see that in 
order to derive the estimate of the scattering amplitude A^, it suffices for us to 
derive the corresponding estimate of ||<l>||rpjj-i/2. To that end, we shall need the 
following lemma, whose proof can be found in [^. 
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Lemma 3.6. Let Br he a central ball of radius R such that Ds (s -Bj?. Then the 
solutions Es,Hs G to ( 1 - 14 ) verify 

/ aiEs ■ rfx + / aoEs ■ E 5 rfx 


' Ds\Ds/2 


>D. 


s /2 


'OBr 

f SR 


{u X EJ) ■ {u X u X H^) da^ + ^ {u x E^) ■ {u x u x H^) 

JdDs 


fz/xE 


'dDs 


■ (u X u X H*) dar- 


dax 

(3.41) 


Remark 3.3. It is remarked that the last two terms in the RHS of (3.41) in 


were 


3? 


'OBr 


{u X E*) ■ (v X ly X Elf )dax + 3R / (z^ x E^) ■ (u x u x H*) dax 


'OBr 


Here, we modify the two terms for the convenience of the present study. 

Dehne a boundary operator A such that 

A(z/ X E^la^J = z/ X Hs\ 9 Br : ^ TR^^^OBr). (3.42) 

It is well-known that A is a bounded operator in (cf. 


using Lemma 3.6, one can show that 


10 


31 


By 


Lemma 3.7. Let E^ and H 5 be solutions to the system (Llf); where Ds is the 

(3.43) 


virtual domain described at the beginning of this section, then we have 


\Esrdx<C6 


i-t 


Di\Ds/2 


\iy xE 


11TH-(aufi) +11 ^ 11 th- 1/2 (az));, 


and 


lEdrdx <C5 


2-t 


'd-s\Ds/2 


lu X Etlp _i/2. 


(3.44) 

where the constant C depends only on D, u and R. 

Proof. First, by (3.41) and the fact that A is bounded in TH^.y^(c}i?H), there holds 

^+||2 


' Ds\Ds^2 


aiEs ■ E^ dx < C*||z/ x + Mi + M2, 


where by using (2.8) and (3.22) one further has 


Ml — 

3? / (z/ X E*) • (z/ X z/ X H+) 

dax 

— 

3R / E* • (z/ X H+) 

dax 


JdDs 

+ 


JdDs 

+ 


< - 
UJ 


ETzi 




(z/x X Vx X Vx X Alsau5i,[a](x)) dax 


+ C6^\\a\ 


TH-U2(az)), 
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and by using (|3.21[) one further has 


M2 — 

3? / (z/ X E+) 

■ [y X V X H*) dax 


[ E+ 

■ (z/ X H*) dax 


JdDs 

+ 


JdDs 

+ 




Combining those with (3.9), (3.11), together with the use of the dehnition of cx/ in 
(3.7), we can complete the proof. □ 


It is remarked that in the estimate (3.44), the dependence on the artihcial R of 
the generic constant C can obviously be absorbed into the dependence on D. 

In what follows, we shall estimate ||‘h||rpjj_i/ 2 (g/) and separately. 

Clearly, it suffices to estimate ||4)||j:^-i/2(5/)3 and ||d)||j:^-i/ 2 ( 5 au 56 ) 3 , respectively. We 
recall that the -norm of the function ‘h(x) is dehned as follows 


II^IIh-i/2(9d) 3 = sup 


'dD 


<l)(x) ■ (y9(x) d(Ti 


Moreover, the H ^/^(S'‘^)^-norm of <h for c G {/, a, b} is given as 


|4>l 


H-1/2(5c)3 : = 


sup 


p/2 


(S=)3 


<1 


$(x) ■(y9(x) dai 


(3.45) 


(3.46) 


where denotes the set of i7^/^(S''^)^-functions which have zero extensions 

to the whole boundary dD. We refer to (^26,35 for more relevant discussions on 
the Sobolev spaces. To proceed, we shall hrst establish the following important 
auxiliary Sobolev extension result, which is a localized version of Lemma 3.4 in [^. 

Lemma 3.8. Suppose D is a simply connected domain with a -smooth boundary 
dD and (s dD is an open surface lying on dD. Let D' he a simply connected 
domain with a -smooth boundary dD' which satisfies ST\dD' = 0 and D' C D. 
Then for any G there exists W G H‘^{Dfi such that 

z/ X W = 0 on S"', 

z/ X (z/ X (V X W)) = z/ X (z/ X -0) on S"^, 

||W||j^2(£,)3 < C||V'||h1/2(5c)3, 

W = 0 in D', 

where C is a constant depending only on D. 

Proof. By zero-extending to dD, we hrst have by Lemma 3.4 in that there 
exists Wo G H‘^{Dfi satisfying 

X Wo = 0 on Sfi 

z/ X (z/ X (V X Wo)) = v X [v X Ip) on S"^, 

I|Wo||h2(£,)3 < C||'0||h1/2(5c)3, 
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where C is a constant depending only on D. Let 6 = supply'll!) denote the snpport 
of the fnnction 'ip. By definition there holds © (s S^. Next, we nse exactly the 
same strategy as that in Addendnm of Theorem 14.1 in [^. We cover © with 
hnitely many, say 91, small enongh regions snch that (U^i H© C 

Then in each ©jH/l there exits an appropriately chosen G 
(see formnla (9) in Addendnm in [^) snch that 

zy X Wj = 0 on n D) n 

z/ X (z/ X (V X Wj)) = u X {v X Ip) on d{Gj fl H) fl S"^, 

6j)nD)3 — ^I|V'IIh1/2(Sc:)3, 

where C* is a constant depending only on D. Finally, by setting W := 

(U^i hlT* and W = 0 in H \ FlD), one can complete the proof. □ 


We proceed with the proof of Theorem 3.1[ 


Lemma 3.9. Let <F he defined in (3.30), where (E^, H^) are the solutions to (1-14) 
with the corresponding e^, ps and as given by (3.1). Then there hold 




(3.47) 


and 


II*^IIth-i/ 2(5 c) < Ch ||E5||£,2 (^c\^^^2 )3, ce{a, 6}, (3.48) 

where (3 = min{l, —1 + r + s, —1 + f + s} and (3' = min{l, —2 + r + s,—2 + t + s}. 

Proof. It snffices to show that the same estimates in ( 3.47| and ( |3.48[ ) hold for 
H-i/ 2 i^scp, c G {f,a,b}. For any test fnnction ip G c G {/, a, 6}, 

we introdnce an anxiliary fnnction W G which satishes the conditions in 

Lemma 3.8 {D' is chosen to be D' = ZI1/2 U {D \ D^)). First, it follows from the 
properties of W that 




<l)(x) • '0(x) dai = - (z/ X E 5 ) • (z/ X (z/ X ^p)) da^ 


'sf 


= - E5 ■ (z/ X (V X W))da^ - 
Jsi 


'dDf\Sf 


E5 ■ (z/ X (V X W)) das; 


>dDf 


E^ ■ (z/ X (V X W)) da^. 


From its constrnction, it is readily seen that u x W|a£)c = 0. Define 

:= (B^)-^(V X W). 


(3.49) 














REGULARIZED FULL AND PARTIAL EM CLOAKS 


21 


By using (3.16), (3.19), (3.12) and integration by parts, one also has 


/ ‘h(x) • '0(x) dax = —6 / E5 ■ (z/ X Curl W ) dax 

' Sf JdDi 


= -6-^ 
=6-^ 
=6-^ 


/ Es ■ [u X CmlW) dax - W ■ [u x CuAW ]dax 

'dDf ^ ' JdDf 


Df 


W-(vxCu^y±-J ■ (^V X Cu^))di 

E^ ■ fv X (V X W)]d±. 


I _ W- V X CurllV)di-5 / 

Df\Dy2 ^ ^ JDf\Dy2 


For X G \ /I5/2 there holds 


Vx X E5 = iu^iHs, Vx X H5 = -iu(^ei + 


By using change of variables and (3.16) in Lemma 3.2, one can derive that 


|B|B ^Vx X E^ = |B|B x + i— ^E^, 


which holds in H \ Di/ 2 - By combining (3.7), (3.15) and (3.49), one can further 
show that in H \ Zli/2, 


V X Curl W = + z(5*+"-)E5. 

CJ 


(3.50) 


Thus by using (3.50) one has 


/ <i>(x) • V’(x) dax 

Jsf 

=a;2(r^+^+" + ir^+*+"-) / _ W-Esdx-df _ E5 ■ (v x (V x W)) di, 

^ JDf\D^/2 dD/\Di/2 ^ ^ 

which in combination with the fact that 

immediately yields that 

I l>(x) • l/>(x) d(T£ <Cd^\\Es\\L^(Df\Du2)^\\'^\\H^iDf\Du2)3 
<C6 ^||E5||^2(£,/\£)^^2)3|| 
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where /3 = min{l, — 1 + r + s, — 1 + t + s}. By using exactly the same strategy as 
above with and (3.20), one can show that for c G {a, b} 


JS‘^ 

=<5-2 


$(x) ■ dax 

f W-(v X ]d5^-6 


[ _ E5- (Vx (Vx W))di, 


JD-\Dy2 ^ 
which together with the fact 

~ ^ ^^‘^\\^s\\l^{d^\Ds/2)^'’ cG{a,6}, 

readily yields that 

[ $ (x) • (x) <i(T:J <C6 ||E5||j;^2(£)c\7)^^2)3i|'^llH2(D=\I?i/2)3> 

where (3' = min{l, —2 + r + s, —2 +1 + s}. The proof can be completed by noting 
the dehnitions (3.45), (3.46) and the fact that 


— ^llV^llHi/2(aD=)U c e {/,a,6}. 


□ 


Proof of Theorem 3J_^ It is straightforward to see from (3.40) that we only need 
to derive the estimates of ||4>||rpjj-i/2c5/^ and ||4>||rpjj-i/2(^g^y To that end, we have 
by using (2.7), (3.13), (3.35) and (3.36) that 


^ ^ ^ C\\u X V X 

— ^^II^IItH-U 2(5/) + C'^^||^|lTH-U2(a£)) (3.51) 

— ^^II‘^IItH-Y2(5/) + C6 (II® ||TH-i/2(g£)) + 1). 


Next, by combining ( |3.51[ ), ( |3.36[ ), ( |3.43[ ), ( |3.44[ ), ( |3.47 ) and ( |3.48[ ), together with 
the use of the facts that both A and A” are bounded, and (3' < fd (see Lemma 


3.9), one can further deduce that 




which in turn yields (noting that (3' — t/2 > 0) 
ll®llTH-i/2(aD) — 
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Thus one has 


1$ 


llTH-U2(a£)) < C5^ 


(3.52) 


It is remarked that in (3.52), the generic constant obviously depends on i?, but as 
was remarked earlier that such a dependence can be absorbed into the dependence 
on D. Next, by using (3.47) and (3.52), there holds 




and thus 




By plugging ( |3.52[ ) and ( |3.53[ ) into ( |3.40[ ), one finally has ( |3.8[ ). 
The proof is complete. 


(3.53) 

□ 


4. Regularized partial-cloaking of EM waves 

In this section, we consider the regularized partial-cloaking of EM waves by 
taking the generating set Tq to be a fiat subset on a plane in M^. In order to ease 
our exposition, we stick our subsequent study to a specific example considered 
in [^[25] for the partial-cloaking of acoustic waves, where Tq is taken to be a square 
on P2 := {x = {xi,X2,X3) E : X3 = 0}. The virtual domain is constructed as 
follows; see Fig. for a schematic illustration. Let n G be the unit normal 
vector to Tq. Let 

Dq := ro(x) X [x - r ■ n, Xr ■ n], x G Tq, 0 < r < g, (4.1) 

where we identify Tq with its parametric representation ro(x). We denote by 
the union of the four side half-cylinders and D'^ the union of the four corner 
quarter-balls in Fig. Let 

D, := D; U DJ U Dl. (4.2) 

Set S'g and 8“^ to denote 

Si := dDq n dDl, si := dD^ n dD^ 

and 

!)/ ;= L)i u D^. 

The upper and lower-surfaces of Dq are respectively denoted by 

Tq := {x-h g ■ n; X G To} and T^ := {x - g ■ n; x G Tq}. 
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Define ^ := rgUr^. We then have dDq = S'°US'g US'g. Similar to onr notations in 
Section we let 5 G M+ denote the asymptotically small regularization parameter 
and be the virtual domain. Clearly, we have 

dDs = 5° U S] U (4.3) 

In what follows, if g = 1, we drop the dependence on q of Dg,Sq, and S^, 
and simply write them as D, S^, S^, and 8“^. In concluding the description of the 
virtual domain for the partial-cloaking construction, we would like to emphasize 
that our subsequent study can be easily extended to a much more general case 
where the generating set Fq can be a bounded subset on a plane with a convex 
and piecewise smooth boundary (in the topology induced from the plane). 



Figure 2. Schematic illustration of the domain Dq for the regular¬ 
ized partial-cloak. 


m 


We introduce a blowup transformation A which maps Ds to D exactly as that 
We stress that in the blow-up transformation takes the following form 


25 


4l(y) = y := 


T 

eseg 


-|- -|- 626^ 
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where y G -D^ and y G -D, and the three Enclidean nnit vectors are as follows 

ei = (1, 0, 0)^, 62 = (0,1, 0)^, 63 = (0, 0,1)^. 

We are now in the position of introducing the lossy layer for our partial-cloaking 
device 

^^(x) = £i(x) := (5^|B|B~\ ns{x) = ni(x) := (5*|B|B"\ 


cr5(x) = cri(x) := (5*|B|B \ for x G -D^ \ -D 5 / 2 , 


(4.4) 


where B(x) ;= Vx^(x) is the Jacobian matrix of the blowup transformation A. 
It is obvious that 


T 

g _ ^3^3 I _ T _ _ j, 


is a constant matrix and 


Bn=-, 


+ eiB; + 6262 in DJ, 


Bei ^ 01, Be2 ^ 02 in 


(4.5) 

(4.6) 


Thus one has 


|B| = 1/5, in Dl 


Furthermore, for x G one can show by direct calculations that (3.10) and 
(3.11) are also valid for the B defined above. 

Lemma 4.1. Let D\ j = 0,1, 2 be defined at the begin ning of this section. Let V, 
W and V, W be defined similarly as those in Lemma 3.2. Then one has 


(z/x X V(x)) ■ W(x)da£ =J ^ / (z/x X V(x)) ■ W(x)da 5 ;, j = 0,l,2. 


'dDi 


'dDi 


We are now in a position to present the main theorem of this section in quanti¬ 
fying our partial-cloaking construction. 


Theorem 4.1. Let Ds be defined in (|4.2[) with its boundary given by Let 

be the pair of solutions to (l.lf) with {1 2; es, Hs, cr, 5 | c /i^, cr^} 

defined in (1.3) and {Ds\Ds/ 2 't£S) given in (f.f). Let A^(x;p,d) be the 

scattering amplitude of Define 


(3j = min{l, - j + r + s, - j + t + s}, j = 0,1, 2, 
and for e G M+ with e -C 1, 

•= {P ^ Up X ii|| < '^d ■= {d G |d ■ n| < e}. (4.7) 

Ifr, s and t are chosen such that (32 — t/2 > 0, then there exists Jq £ 1^+ such that 
when 6 < 5o, 

||AL(x;p,d)|| + pGSp, d G x G (4.8) 

where C is a positive constant depending on oj and D, but independent of ea, fVa, 
(Ttt and X, p, d. 
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Remark 4.1. Similar to Remark 3.1, one readily has from Theorems l.l| and 4.1 that 
the push-forwarded structure in (1.2) produces an approximate partial-cloaking 
device which is capable of nearly cloaking an arbitrary EM content. The highest 
accuracy of such a construction can be achieved, say for example, by taking r = 
t = 0 and s = 5/2. 


Remark 4.2. Since p ■ d = 0 (cf. (1.7)), one has 

(d ■ n)p = d X (p X n). 


(4.9) 


from which one can infer that the definition of in (4.7) is redundant. However, 
we specihed it for clarity and definiteness. 

4.1. Auxiliary lemmas and asymptotic expansions. We shall follow a similar 
strategy of the proof for Theorem 3T in proving Theorem 4T In what follows, we 
adopt similar notations as those in Section If we let z denote the space variable 
on To, then for any y G dDq, we dehne Zy to be the projection of y onto Tq. We 
also define a(y) := a(y) and y := A(y) G H for y G D^. 

The following lemma is a counterpart to Lemma 3.4 in Section 


Lemma 4.2. Let Ds be described in (4.2) and (4.3). Define 


Algo [a] (x) := p.v. z/x X Vz- X / - Zy)a(y) day, 


(4.10) 


/go 


and 

Define 

and 


Alg 2 [a](x) := p.v. L-x X Vx X / G^(x - y)a(y) day. 


'S2 


>^‘Dja](x) := z/x X Vx X Vx X A‘/,Ja](x), 

/:go[a](x) := z/x X Vzj X Vzj X / - z^)a(y) day. 

Jso 


Then one has 

Jl ) \ Ai;^o[h](k) + Alg2[d](k)+ C>(5||h||TH-i/2(az5)), xgR|, 


Algo [a] (x) -|- d((5||a||rpjj_i/2('gg,^), x G 5° U Sj, 


and 
^Dja](x) = 


(4.11) 


_ / '^so[a](x)-|-(!l(5||a||rpjj_i/2(-gg,^^ 


gr°uri. 


7z/x X Vx X Vx X Ag 2 [a](x) + 0(||a||TH-i/2(aD))> ^ G Sj. 


Proof. By following a similar argument to that in the proof of Lemma 3.4 
compute for x G S'® U S'] that 


(4.12) 
one can 
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and 

(x - y, = (z, - z„ + 0{5), + 0{5). 

One also notes that for x G S'! 

||x - y|| = 5||x - y||, = 1 + 0{5). 

By nsing the above facts, the proof can then be completed by nsing the similar 
expansion method as that in the proof of Lemma |3.4[ □ 

Next, by straightforward calculations, one can show that for described in 


(4.2) and (4.3), there holds the following far field expansion for x G \ D 5 , 


5^Ja](x) = / G^(x - y)a(y)d(Ty = / G^(x - Zy)a(y)dcrji 
JdDs Jso 

+ 6 G^(x - Zy)a(y)dcry + 5 / (y - z^) ■ VG^(x - Zy)a(y)d(Tj; 


'51 


Iso 

l^llTH-U 2 (a_D))- 

With the above expansion, we can show the following theorem 


(4.13) 


Theorem 4.2. Let be the solution to (1-14) with Ds described in (4.2) and 

+ 


(4.3). Define <h(y) := *h(y) = Py x E 5 (y) , then there holds for x G R^\D 

dDs 


(E,-E')(x)= / G',(x-z^)M[nx E'(z)](z^)da, 
Jro 


+0 Mix 


1-1 


I*^IItH-i/2(50) + <^(||‘^||TH-U2 (a£)) + 1) j j , 


where 


with defined by 


-(- 4 +-^ro 


-1 


-^ro ( 4 + -^ro 


-1 


(4.14) 


(4.15) 


:= n X Vzi X / Ga;(zx - Zy)Q{zy) da^ 
Jto 


(4.16) 


Proof. We hrst recall that the solution to (1.14) can be represented by (2.7) and 


(2.8), where the density function a satisfies (2.9). Next, by using the fact that 
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is invertible (cf. |Mp^ ), and using (2.9), (4.11), along with the use of the expansion 
of E* in z, one has by direct calculations that 


^(y) “ 2 


V X E*(z) + 4) 

= (-i + V<S.)"LxE‘(z)] (y) 


(y) + ^('^(ll‘^llTH-U 2 (a£)) + 1 )) 


+ O ( ||‘h||rpjj-l/2(g0) + ^(||‘h||rpjj-l/2(g£,) + 1) 


(4.17) 


for all y e 5'° U S^. Noting that for y G F^, (y — 2n) G F^, we dehne i/’(y) := 
a(y — 2n) for y G F^ and '(/'(y) := ^(y + 2n) for y G F^. By using the fact that 

^y- 2 ny = -^y = “H for y G F\ 

and the dehnition (4.10), one obtains (assuming for a while that x G S'°) 

n X E*(zi) (y) 


a(y - 2n) = V’(y) = (-+M 


so 


+ ^ ( ll^ll ( + <^(||4>||Tjj-l/2(-g^) + 1) 


(4.18) 


for all y G F^. By inserting (4.17) and (4.18) into (4.13), one then has that for 
X G \ )D, 


^Ds[^]i^) = - z^)a(y) day + O (^||x||- (^5(||<l)||TH-i/2(9i5) + 1) 

= ^.(x - z^)(a(y) + Vi(y)) day + +0 (||x||-^ (^(ll^llTH-i/2(aD) + 1))) 

= [ G^(x - Zjj)M[n X E*(z)](y) 

2 ri 

+ ^ (^ll^ll (^II^IItH-Y2(50) + <^(||^||TH-Y2(g^) + 1)^) ) 
where M is dehned by 




foo 


(4.19) 


Since the function 0(z) := n x E*(z) depends only on z G Fq, and hence by 
dehnition (4.16) one readily verihes that 

-^ro[0](z£) = ^>^so[0(y)](z£+ n), iG Fi, 


for 0(y) = n X E*(zy), y G S^. Therefore, we can replace the operator M in (4.19) 


to be 


-1 


:= (-T+-Mrol 


Lo 


‘Lo 


-1 
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and the proof is complete. 


□ 


By using Theorem |4.2 , we next consider a particular case by assuming that 
$ = 0. Physically speaking, this corresponds to the case that is a so-called 
perfectly electric conductor. The result in the next theorem already partly reveals 
the partial-cloaking effect. 


Theorem 4.3. Let Ds be as described in (4.2) and (4.3). Consider the following 
scattering problem 


V X Ef - iuW^ = 0 in \ 

V X = 0 in \ 

z/ X E||_i_ = 0 on dDs, 

subject to the Silver-Miiller radiation condition: 


lim ||x|| 

|x|^cx: 


(H? 


H*) X X - (Ef - E*)) = 0. 


(4.20) 


(4.21) 


Let E^(x;E*) be the corresponding scattering amplitude to (4.20)-(4.21). Then 
there holds 


E^(x;E*) + - 


E Y{-{k)Yp-{Zy)\Zy\ 


To 


[nxE\z)]{zy) da^. < (4.22) 


where M is defined in (4.15) and C depends only on uj and D. Furthermore, if 
there holds 

||n X pII < e, e 1, 
then for sufficient small 6 G M+, one has 


|EL,(*;E' 


< C(e + i), 

where C depends only on oj and D. 

Proof. We start with the following addition formula (see, e.g., |M]) 


(4.23) 

(4.24) 


1 




2n + 1 






,fn 


n+1 ' 


9 


(4.25) 


" " n=0 m=—n 

where (g, i?) and (g', -d') are the spherical coordinates of x and y, respectively; 

and Yffi is the spherical harmonic function of degree n and order m. For simplicity, 
the parameters (g,^,!?) and shall be replaced by (||x||,x) and (||y||,y), 

respectively. It then follows by (4.25) that 


471 


|x-y|| = ||x|| 1-- E U’‘(*)L"‘(y) 


+ 0(||x 


1 - 1 ) 


(4.26) 


m=—l 


The solution E^ in (4.20) and (4.21) can be represented by 


E' = E‘ + ylS,[a 


m 


\Ds, 
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with a satisfying 




[a](y) = -I^y X E*(y), yedDs. 


Then one can easily derive from the expansion formula (4.14) with <l> = 0 that 

(E^ - E*)(x) = 2 [ - Z 5 )M[n x E'(z)](z^) da,. + 0(||x||-i5), x G \ Ds, 

Jfo 

(4.27) 

where M is dehned by (4.15). By combining (4.27) and (4.26), we then have 

1 


E^(i;E-) = -- 


-iuJ^ E Xr(x)Yi’"(zy)|zjj| 


To 


[n X E*(z)](zy) da,- + 0{6), 

(4.28) 


which readily proves (4.22). Next by the dehnition of M in (4.15), one can show 
that 

||M|nxE-(z)l|| TH-i/2(ro) < C||nx 

E*(z)||TH-U2(ro)5 

where C depends only on Tq and u. This together with (4.28) further implies that 
||Ej„(x;E‘)|| < C(||n X E'(z)||TH-./>(r.) + i)- 


(4.29) 


Finally, by inserting the condition (4.23) into (4.29) and direct calculations, one 

□ 


can easily show (4.24). 


The proof is complete. 

4.2. Proof of Theorem |4.1| , In this section, we present the proof of Theorem |4.1[ 
which follows a similar spirit to those of Theorems 3T and |4.3[ The major idea 
is to control the norm ||‘h||rpjj-i/ 2 ('g^), which was taken to be identically zero in 
Theorem 14.31 

Lemma 4.3. Let (E^, H^) be the pair of solutions to h with jn- es, las, as} C 
£ 5 ,/X 5 , cr^} defined in (1.3) and {Ds\Ds/ 2 ]£ 5 , t^s, 0 's} given in (ffi). Then there 
holds the following estimates for j = 0,1, 2, 


'DI\Ds/2 


lE^II dx <C5^ ( ||z/ X E^ + ||^|ItH-iG( 50 ) + <^||^|lTH-iG(a£)) j, 


where the constant C depends only on R and uj. 


the present lemma. Then one has 


Proof. We hrst note that (|3.41|) still holds for the scattering problem described in 


/ _ (TiE^ ■ E^dx < C||i/X E7||.^jj- 

'Ds\Ds/2 
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where by using (2.8) and (4.12) one further has 


M.]^ — 

n ! (i/ X E‘) ■ (i/ X z/ X H+) 

dax 


^ [ E/ ■ (z/ X H+) 

dax 


JdDs 

+ 


JdDs 

+ 


1 

< - 
OJ 


Iso 


E*(zx) ■ -^so[a](x) dax + (!2(5||a||rpjj-i/2^gj)^), 


and by using (|4.11[) one further has 


M 2 — 

uj (i. X E+) 

■ {y X V X H*) dax 

— 

3fJ / E+ 

■ (z/ X H*) dax 


JdDi 

+ 


JdDs 

+ 


< 5 


Iso 


Al5o[a](x) ■Yl\7,x)dax + C>((5^||a||.pH-i/2(g^0. 

Combining the above estimates with (4.5), (4.6) and the dehnition of ai in (4.4), 
we can complete the proof. □ 

Proof of Theorem 4^. First, by using (4.14) and (4.7) one obtains 

P) d)ll ^ + ||$||th-U 2 ( 50 ) + <^(||^||TH-U2(a£)) + 1)^, (4.30) 

where *h(x) = ‘h(x) := p x (x) for x G dDg. The following estimate can be 

+ _ 

obtained by using the result in Lemma 4T and a completely similar argument as 
that in the proof of Lemma 3.9 

II^IIth-i/2(5j) < j = 0,l,2, (4.31) 

where {3j = min{l, —j + r + s, —j + t + s},j = 0,1, 2. On the other hand, by using 
(4.17), Lemma 4.3| and (4.30), one has 

[ + ( 4 . 32 ) 

+ II®IItH-iG( 50 ) + <^||‘^||TH-U 2 (g£)) + e + , j = 0,1,2. 

Inserting (4.32) back into (4.31), one can show 

II ^ IITH- 1 / 2 ( 52 ) <05^^ ^ ^ (^II^IItH-i/ 2(S0) + ^ ll^llTH-i/2(aT>) 

+ 11411 ™-./^.,+«“'' +' 5 '''). 1 = 0 . 1 . 2 . ( 4 . 33 ) 

Noting that (32 < jdi < /3o and /?2 — ^/2 > 0, one has for 5 G M+ sufficiently small 
that 

Inserting (4.34) back into (4.33) (for j = 0), there holds 

+ II^IIth-i/2(so) + (4.35) 


|‘h||TH-i/2(50) — (||®IItH-i/2(S0) 
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By /32 — ^/2 > 0, it is directly verified that Po — t/2>l and hence 

2(/3o-t/2)-l>/?o-t/2. (4.36) 

Using (4.35) and (4.36), one readily infers for 5 G M+ snfficiently small that 

-1/2(50) <C'(<5^°-*/2 + e). 


|4>l 


TH- 


(4.37) 


Now by inserting (4.37) into (4.34), one has 

< C'5/32-t/2-l/2^^1/2 ^ ^ 1 / 2 ) < + 52(/32-t/2)-l ^ ^4 33 ) 


Finally by plngging (4.37) and ( |4.38| ) into (4.30), we arrive at ( |4.8[ ). 
The proof is complete. 


□ 
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